RELATIVISTIC ENTROPY AND RELATED BOLTZMANN KINETICS 
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It is well known that the particular form of the two-particle correlation function, in the collisional 
integral of the classical Boltzmman equation, fix univocally the entropy of the system, which turn 
out to be the Boltzmann-Gibbs-Shannon entropy. 

In the ordinary relativistic Boltzmann equation, some standard generalizations, with respect its 
classical version, imposed by the special relativity, are customarily performed. The only ingredient 
of the equation, which tacitly remains in its original classical form, is the two-particle correlation 
function, and this fact imposes that also the relativistic kinetics is governed by the Boltzmann- 
Gibbs-Shannon entropy. Indeed the ordinary relativistic Boltzmann equation admits as stationary 
stable distribution, the exponential Juttner distribution. 

Here, we show that the special relativity laws and the maximum entropy principle, suggest a 
relativistic generalization also of the two-particle correlation function and then of the entropy. The 
so obtained, fully relativistic Boltzmann equation, obeys the H-theorem and predicts a stationary 
stable distribution, presenting power-law tails in the high energy region. The ensued relativistic 
kinetic theory preserves the main features of the classical kinetics, which recovers in the c — > oo 
limit. 



I. INTRODUCTION 

In experimental high energy physics, the power law 
tailed probability distribution functions have been 
observed systematically (plasmas [1| , cosmic rays 0, |3l , 
particle production processes d, @ etc). 

A mechanism, frequently used to explain the occur- 
rence of non exponential distributions, is based on cer- 
tain non-linear evolution equations, rnainl y co nsidered in 
the Fokker-Planck picture 0, i, i, El Q IH, 

[l6l.[l7l. isl-irpl. but recently also in the Boltzmann picture 



ToT 20j . Clearly, the correctness of the analytic 
expression of a given distribution, used to describe a sta- 
tistical system, is strongly related to the validity of its 
generating mechanism. 

The classical Boltzmann equation, due to the partic- 
ular form of the two-particle correlation function, in the 
collisional integral, fix univocally the entropy of the sys- 
tem, which turn out to be Boltzmann-Gibbs-Shannon en- 
tropy. The latter entropy imposes the exponential form 
to the probability distribution function, emerging as the 
stationary and stable solution of the equation. 

In the ordinary relativistic Boltzmann equation [2l| . 
some standard generalizations, with respect its classical 
version, imposed by the special relativity, are customarily 
performed. The only ingredient of the equation, which 
tacitly remains in its original classical form, is the two- 
particle correlation function, and this fact imposes that 
also the relativistic kinetics is governed by the classical 
Boltzmann-Gibbs-Shannon entropy. Indeed the ordinary 
relativistic Boltzmann equation admits as stationary sta- 
ble distribution, the exponential Juttner distribution. 
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In order to obtain a relativistic Boltzmann equation 
admiting as stationary and stable solution a probability 
distribution, different from the exponential one, the only 
possibility we have, is to modify the expression of the 
classical two-particle correlation function. As a conse- 
quence a modification of the system entropy emerges, and 
this important fact suggests that such modification must 
be confined possibly within the ordinary physics (spe- 
cial relativity, maximum entropy principle) and without 
invoking any additional assumption or postulate. 

Main goal of the present paper is to show that it is 
possible to obtain a fully relativistic Boltzmann equation 
admitting as stationary and stable solution a distribution 
function presenting power-law tails. 

The new relativistic equation, can be obtained start- 
ing from the ordinary relativistic Boltzmann equation 
[2ll | , by properly generalizing the two-particle correlation 
function in the collisional integral. 

Such generalization, can be obtained in a self- 
consistent manner, by employing (i) the special relativity 
laws and (ii) the maximum entropy prin ciple, and with- 
out invoking any extra principle [1^, [ij, [H [l^ . 

The new relativistic entropy and relevant probability 
distribution function have very simple expressions, which 
result to be one-parameter (light speed in dimension- 
less form) generalizations of the corresponding classical 
concepts, just as happen for all the physical quantities 
(energy, momentum etc) in special relativity. Clearly, 
the relativistic entropy and the relativistic distribution 
function, in the classical limit, reduce to the Boltzmann- 
Gibbs-Shannon entropy and to the Maxwell-Boltzmann 
distribution respectively. 

Within the present theoretical framework, the power- 
law tails in the distribution functions, in high energy 
physics, emerge as a purely relativistic effect (those tails 
in the classical limit c — s- cx) deform into exponential 
tails). In other words the power-law tails represent the 
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signature of the relativistic nature of the system. 

The paper is organized as foUows: 

In Sect. II we consider briefly the basic concepts of 
relativistic dynamics by using dimensionless variables. 

In Sect. Ill we consider the K-differential calculus by 
giving particular emphasis to its physical origin. 

In Sect. IV we present the main properties of the func- 
tion K-exponential. 

In Sect. V we present the main properties of the func- 
tion K-logarithm. 

In Sect. VI we introduce the velocity representation of 
the functions K-exponential and K-logarithm. 

In Sect. VII we consider Boltzmann relativistic kinet- 
ics and discuss in detail the relativistic generalization of 
the classical Boltzmann-Gibbs-Shannon entropy. 

Finally in Sect. VIII we report some concluding re- 
marks. 



II. RELATIVISTIC DYNAMICS WITH 
DIMENSIONLESS VARIABLES 

In the present section, for simplicity of the exposition, 
we consider briefly, the relativistic dynamics in one spa- 
tial dimension. Starting from the velocity v, the momen- 
tum p and the total energy E of a particle of rest mass 
m we introduce the dimensionless velocity u, momentum 
q, and total energy £, through 



V 

u 



P 

mq 



= KC = IK 



(2.1) 



being c the light speed and < c a reference velocity. 
The condition —c<v<c implies that —1/k < u < 
so that 1/k represents the dimensionless light speed. Al- 
ternatively 1/k can be viewed as the refractive index of a 
medium in which the light speed is v^, . From the disper- 
sion relation E'^ = m^c^ +%P'c^ follows that E > mc^ and 
consequently it holds S > 1/k^. Clearly the condition 
—00 < p < +O0 implies — oo < q < +oo. For a particle 
at rest it results E{0)/£{0) ~ mc^K^ and after posing 
E{0) = mc^ we obtain £{Q) — 1/k^. Then we can con- 
clude that 1/k^ represents the dimensionless rest energy 
of the particle. It is important to note that the classical 
limit limc^oo in dimensionless variable transforms into 
the limit lim„^o- 

In relativistic dynamics any physical quantity or vari- 
able /i, is linked through an invertible function with the 
particle momentum i.e h — h(q), q = q{h). The velocity 
u and the total energy £ are given by 



£{q) ^ \^l + nW , 



(2.2) 
(2.3) 



while the kinetic energy >V = f — 1/k^>0 assumes the 
form 



w(g) = ^Vi + «V-^ 



(2.4) 



Clearly we can introduce further function h = h{q) 
which define new variables. For instance the Minkovski 
rapidity ■(/;, is defined through [3 — tanh'0 with (3 — v/c. 
After introducing the new rapidity variable p through 
"0 = up, we can express it in terms of dimensionless mo- 
mentum as follows: 
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p{q) = — arcsinh(Kg) 



(2.5) 



The above definition of the rapidity variable is justified 
by the fact that in the classical limit it results: p = q = 
u ^ V2W. 

Let us consider in the inertial one-dimension spatial 
frame E, two identical, non interacting, and free particles 
A and B, of rest mass m. We suppose that the particle 
A moves toward right while the particle B moves toward 
left. The dimensionless variables for the particle A are 
given by g^, Va, £a, VVa and Pa- The absolute values of 
the momentum, velocity and rapidity of the particle B 
are given by , Vb, Pb, while its total energy and kinetic 
energy are given by £b, Wb respectively. 

In the rest frame E' = of the particle B, which 
moves toward left, with velocity Vb with respect E, the 
above considered dimensionless variables, for the particle 
B are given by q'^ —Q,v'g— 0, £'g — 1/k^, W'b ~ and 
p', =0. 

According to the Lorentz transformations, the dy- 



namic variables of the particle A i.e. q'^ 



£', m 



and p'a, in the frame E', are direct compositions of the 
corresponding variables of the particles A and B in the 
frame E. For instance the momentum composition law 



9a = Qa ( 



(2.6) 



with 



Qa® Qb = Qa y/l + K^ + qB y/l + K^ 



(2.7) 

follows directly from the Lorentz transformation of the 
four-vector energy-momentum. 

Regarding the composition law of an arbitrary variable 
h{q), after posing hj^ = h{qA), = h{qB) and /i^ = 
h{q'ji^, one obtains 



with 



h'^ = hA®h^ 



h^®h^ =h{q^®q^) 



(2.8) 



(2.9) 

so that the composition law of any variable can be de- 
duced from the momentum composition law. 

In particular starting from the general composition law 
defined in Eq. (|2.9p and after taking into account Eqs. 
()2.2|) - ()2.5p we can obtain the composition laws for the ve- 
locity, total energy, kinetic energy and rapidity as follows 



~ Ua ® Ub , 
£a ^ ^A ® £b ; 
Wa^Wa®Wb , 

Pa^ Pa® Pb , 



(2.10) 
(2.11) 
(2.12) 
(2.13) 
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being 



11a VP Ub — 



Ua + Ub 
1 + K^UaUe 



(2.14) 



1 



£a® £b^ k'EaEb + ^^{n^£l~l){K^£l-l), (2.15) 

W^d Wb = Wa + Wb + h^WaWb + 

+ y^>V4ni(2 + K2>V^)(2 + K2>v^) ^ (2.16) 

Pa® Pb ^ Pa+ Pb ■ (2-17) 

It is remarkable that the momentum composition law 

9a ffi fe, is a generalized sum having the properties: i) 
it is associative, ii) it is commutative, iii) the is its 
neutral element, iv) the opposite element of q is —q. Also 

the velocity composition law Ua ® Ub, results to be a 

K 

generalized sum, while the rapidity composition law Pa © 
Pb is the ordinary sum. 

Regarding the total energy composition law £a ffi £b 
we note that i) it is associative, ii) it is commutative, 
iii) it has as neutral element the quantity 1/k^, iv) it 
does not exist the opposite element of £. Due to the 
latter property, the total energy composition law is not 
a generalized sum. Also the kinetic energy composition 

law Wa ® yVs is not a generalized sum because it does 
not exist the opposite element of W. 

III. THE K-DIFFERENTIAL CALCULUS 

A. Lorentz invariant integration 

Let us consider the following integral in the three di- 
mension momentum space, within the classical physics 
framework 



(3.1) 



being F an arbitrary function and = men. In dimcn- 
sionless variables, the latter integral becomes 



Xo = / d\F 



(3.2) 



Whenever F depends only on q = |q|, the above integral 
can be reduced to the following one dimension integral 



Jo = / dq 47r q^ F 



(3.3) 



In the framework of a relativistic theory it is well 
known that the three-dimension integral (j3.1[) must be 
replaced by the four-dimension Lorentz invariant integral 



mc 



2e{pQ) 5{p^'p^-m^^) F . (3.4) 



After introducing the dimensionless variables according 
to qf^ = p^/p* = (gojq) = (Kf,q), the latter integral 
becomes 

In= Jd\2e («:go) S {t^q^'q^, ^ I) F , (3.5) 

or alternatively 

j d^q J d{K^8) 2e{K^£) S {k^8^ - n^q^ - l) F 

(3.6) 

with q = |q|. 

In order to reduce the integral (|3.6p . in a three- 
dimension integral, we introduce the new integration 
variable Z = k'^E'^. After observing that d{K^£) ~ 
dZ/2^fZ and e{K'^£) = 9{Z) one obtains 



1,^1 d\ £^ ^2S{Z- -I) F (3.7) 



and finally 



dh^ 



F . 



(3.8) 



The latter three-dimension integral in dimensional vari- 
ables becomes 

^[fj^^F (3.9) 

J Pi Po 

We remark that in Eq. p.4p the integration element 
d'^p, is a scalar because the Jacobian of the Lorentz trans- 
formation is equal to unity. Then Irei transforms as 
F. For this reason in Eq. p.Sp the integration element 
d'^q/ \/l + n^q^ is a scalar. 

Whenever F depends only on g = |q|, the integral (|3.8p 
can be reduced to the following one dimension integral 



dq 



Airq^F , 



10 y^l + n'^q^ 
which after introducing the K-differential 

dq 



d^q = 



assumes the form 



d^q An q^ F 



(3.10) 



(3.11) 



(3.12) 



In the classical limit k ^ 0, the latter integral reproduces 
the corresponding classical one given by Eq. p.3p . 

We focus now our attention to the classical and rela- 
tivistic expression of the one dimension integrals given by 
p.3p and (|3.12p respectively. One immediately observes 
that the relativistic integral is obtained directly from the 
classical one, by making the substitution dq — > d^f?. 
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The equivalence of the Lorentz invariant integration 
in four-dimension J d^p, with the ^-integration in one- 
dimension Jp°° dt^q, according to 

/"OO 

d^p e (po) S {p^p^ ~ Tr?^) F{p) ^ d^q q^F{q) , 

Jo 

(3.13) 

permits to explain better the relativistic origin of the k- 
integration. 



B. Kinetic Energy and Work 

In classical mechanics the kinetic energy W is defined 
as the work produced by the external force f 



Jo 



(3.14) 



After taking into account the Newton law and using di- 
mensionless variables the above definition assumes the 
form 



>Vo(<z) ^ qdq 



(3.15) 



and yields Woiq) = q'^/2. 

In special relativity, by using the same definition, the 
kinetic energy is given by 



Jo Jo Jo 

V dp — ' 







v^ + pV^ 



vdp 

= dp . (3.16) 



In dimensionless variables, the kinetic energy becomes 

(3.17) 



>V.(.)-/%^=i^ ' 
Jo VI + K^q^ 

and after introducing the K-differential can be written in 
the form 



W«(9) ^ q d^q 
Jo 



(3.18) 



obtaining in this way the well known relativistic expres- 
sion W^{q) = (Vl + K2g2_ 

It is remarkable that the replacement of the ordinary 
integration by the K-integration, in the classical definition 
p.lSp . is sufficient to recover the relativistic expression 
of the kinetic energy. 

Eq. p.lSp linking W and q through the K-integral, can 
be written also in the following differential form 



Eq. (|3.19[) after integration with the condition Wk,{0) = 
0, yields the relativistic expression of the kinetic energy. 
In the K ^ Q limit the differential equation p.l9p . re- 
duces to the classical one 



dq 



(3.21) 



C. Momentum generalized sum and K-differential 
calculus 

Let us consider in the inertial one-dimension spatial 
frame S, two identical, non interacting, free particles A 
and B of rest mass m which move toward right. We sup- 
pose that the dimensionless momenta for the two particle 
are given by = 9 + dq, qg = q respectively. 

In the rest frame S' = , of the particle B the particle 

K 

momenta are given by = and q'^ ~ {q + dq) © (—'3') 
respectively. 

The infinitesimal difference of the particle momenta, 
in the frame S is dq, while in the frame E' is given by 



d^q = (q + dq) {-q) 

K, 

= {q + dq) Qq , 



and results to be 



d^q 



dq 



q^ 



(3.22) 
(3.23) 



(3.24) 



In order to better understand the origin of the ex- 
pression of the K-differential, we recall that the vari- 
able g is a dimensionless momentum. Then the quan- 
tity 7(g) = -\- q^ is the Lorentz factor in the mo- 
mentum representation, so that the K-diffcrential can be 
written as 



d>^q 



dq 
7(9) 



(3.25) 



In other words if dq is the infinitesimal difference of two 
particle momenta in the inertial frame S, this difference 
if observed in rest frame S' of one of the two particles, 
becomes dK,q and results to be contracted by the Lorentz 
factor. 

The K-derivative of the function /(<?), is defined 
through 



^^limM. 



q 



zQ q 



f{q) _ f{q + dq)-f{q) 
{q + dq)eq 



(3.26) 



dfcq 



involving the K-derivative 

d 
dttq 



d 

dq 



(3.19) 



(3.20) 



We observe that df{q)/d^q, which reduces to df{q)/dq 
as the deformation parameter k — > 0, can be written in 
the form 



dfjq) 

dK q 



K? q^ 



dfjq) 
dq 



(3.27) 
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From the latter equation follows that the generalized 
derivative obeys the Leibniz's rules of the ordinary 
derivative. After introducing the Lorentz factor ^(q), the 
derivative operator can be written also in the form: 



d r \ d 

-1 = 7(9) -T- ■ 

d^q dq 
The K-integral is defined through 



d^q f{q) = 



dq 
7(g) 



(3.28) 



(3.29) 



and obeys the same rules of the ordinary integral, which 
recovers when k ^ 0. 



IV. THE K-EXPONENTIAL FUNCTION 
A. Definition 

In the present section, the independent variable is a 
dimensionless momentum and is indicated hy x, y or 
z. We recall that the ordinary exponential f{x) — 
exp(a;) emerges as solution both of the functional equa- 
tion f{x -\-y) = f{x)f{y) and of the differential equation 
{dldx)f{x) = fix). 

The question to determine the solutions of the gener- 
alized equations 



fix ®y) = fix)fiy) 



dfix) 
dKX 



fix) , 



(4.1) 



(4.2) 



reducing in the k — > limit, to the ordinary exponen- 
tial, naturally arises. The latter two equations admit the 
same solution, which represents an one-parameter gener- 
alization of the ordinary exponential. 

Solution of Eq. CT71); We write this equation explicitly 



/ (xVi + n^y^ + yVi + ^^x^ ) = fix)fiy) , (4.3) 

and after performing the change of variables fix) = 
exp((7(Kx)), zi = Kx, Z2 = Ky it transforms into the 
following equation 

9 (^^i\/l + 2| + 22^1 + 2? j =g(zi)+5(z2) , (4.4) 

which admits the solution g(x) — Aarcsinhx. Then, 
it results that fix) = exp(A arcsinh kcc). The arbi- 
trary constant A can be fixed through the condition 
Imii^^Q f ix) — exp(a;), obtaining A = l/n. Therefore 
fix) assumes the form fix) = exp^(a::) being 



exp^ ix) — cxp ( — arcsinh kx 



(4.5) 



Solution of Eq. ^7^: After performing the change of 
variable p = arcsinh ko;, it obtains d^a; — dp and Eq. 
(14.911 assumes the form 



dl 
dp 



= f ■ 



(4.6) 



From the solution of the latter equation with the condi- 
tion /(O) = 1, follows immediately that fix) — exp^(a:) 
with 



(4.7) 



exp^ (x) = exp — arcsinh kx 



By taking into account that arcsinh a; — ln(\/l + x"^ + x) 
we can write exp^ (x) in the form 



exp, 



. (x) = (^\/ 1 + K^X^ + KX^ 



1/k 



(4.8) 



which will be used in the following. We remark that 
exp^(a;) given by Eq. (|4.8p . is solution both of the Eqs. 
(14. 1|) and (14. 2p and therefore represents a generalization 
of the ordinary exponential. 

In particular according to Eq. (|4.2p the ^-exponential 
is defined as eigenfunction of the K-derivative i.e. 



d 



exp„(x) = exp„(a;) 



(4.9) 



B. Basic Properties 

From the definition (|4.8p of exp,,(.-E), follows that 

expg(2;) = limexp^(x) = exp(x) , (4.10) 
exp„^(j;) = exp^(a;) . (4.11) 
Like the ordinary exponential, exp^ (x) has the properties 

(4.12) 
(4.13) 
(4.14) 



exp,(x) e C°°(R), 
d 

— exp,,(a:;) > 0, 
dx 

exp^i-00) = 0+, 
exp«(0) = 1, 
exp«(+oo) = +00, 
exp^ix)exp^i-x) = 1 



(4.15) 
(4.16) 
(4.17) 



The property (|4.17p emerges as particular case of the 
more general one 

exp^(a;) exp^iy) = exp^(x ® y) . (4.18) 

Furthermore exp^(a::) has the property 

[exp^(x)]'' = exp^/^(ra;) , (4.19) 

with r G R, which in the limit k ^ reproduces one well 
known property of the ordinary exponential. 
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We remark the following convexity property 



2 exp^(x) > ; X G R , 



(4.20) 



holding when < 1. 

Undoubtedly one of the more interesting properties of 
exp^(x), is its power law asymptotic behavior 



exp^(a;) ~ 2kx 

X — >±oo 



±1/|k| 



(4.21) 



It is remarkable that the first three terms in the Tay- 
lor expansion of exp^ (a;) are the same as those of the 
ordinary exponential 

exp,(a:) = l + x+y + (1-At2)|! + ... . (4.22) 

This latter result is a particular case of a more general 
property of the relativistic dynamics. Indeed, the Taylor 
expansion up to the second order, of any relativistic for- 
mula, coincides with the corresponding classical formula. 



V. THE K-LOGARITHM FUNCTION 
A. Definition 

The function hii^(x) is defined as the inverse function 
of exp^ (x) , namely 

lnK(exp^ x) = exp^(lnK x) = x , (5.1) 

and is given by 



hii^{x) — — sinh(Klnx) , 

K 



or more properly 



\n^{x) 



2k 



(5.2) 



(5.3) 



B. Basic properties 

It results that 

Ino(x) = lim \i1k{x) = ln(x) , (5-4) 

(5.5) 



ln_K(x) = liif,{x) 



The function ln„(x), just as the ordinary logarithm, 
has the properties 



Injx) eC°°(R+), 


(5.6) 


ln^{x) > 0, 

ax 


(5.7) 


lnJO+) = -oo, 


(5.8) 


ln«(l) = 0, 


(5.9) 


ln„(-|-oo) = -l-oo. 


(5.10) 


lnK(l/x) = -lnK(x) . 


(5.11) 



Furthermore InK(x) has the two properties 

ln,(x'^) =rln,,(x) , (5.12) 
ln«(x?;) = ln«(x) ® ln«(y) , (5.13) 



with r G R. Note that the property (|5.1ip follows as 
particular case of the property (|5.12[) . 

We remark the following concavity properties 



X 

-7-^ X \ii^{x) < 
ax^ 



-—r \n^{x) < , 
dx'^ 



(5.14) 
(5.15) 



A very interesting property of this function is its power 
law asymptotic behavior 



ln«:(x) 
ln«;(x) 



1 



x^o+ 2|k| 
1 



■>+oo 2 k 



(5.16) 
(5.17) 



After recalling the integral representation of the ordi- 
nary logarithm 



1 1 
ln(x) = - / -dt 



(5.18) 



one can verify that the latter relationship can be gen- 
eralized easily in order to obtain ln„(x), by replacing 
the integrand function ya{t) = by the new function 

y^it) =1^^^*^, namely 



ln«(x) 



1 



1 



2 Ji/x 



dt . 



(5.19) 



The first terms of the Taylor expansion related to k- 
logarithm, are 



ln.(l + x)=.T--+(^l + -l--... . (5.20) 



C. Tlie lnK,{x) as solution of functional equations 

First functional equation: The logarithm f{x) — ln(x) 
is the only existing function, except for a multiplicative 
constant, which results to be solution of the function 
equation f{xiX2) = /(xi) -I- /(X2). Let us consider now 
the generalization of this equation, obtained by substi- 
tuting the ordinary sum by the momentum generalized 
sum 



fixiX2) = /(xi) ® /(xa) 



(5.21) 



We proceed by solving this equation, which assumes the 
explicit form 



/(X1X2) - /(Xi) Vl + «2/(x2)2 
+ /(X2)V1 + «V(X1)' 



(5.22) 
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After performing the substitution f{x) = sinh Kg {x) 
we obtain that the auxiUary function g{x) obeys the 
equation g{xiX2) = g{xi) + g{x2), and then is given 
by 9{x) — Ahix. The unknown function becomes 
f{x) — sinh(Klna;) where we have set A = 1 in or- 
der to recover, in the Umit k ^ 0, the classical solution 
f{x) = ln(a;). Then we can conclude that the solution of 
Eq. (|5.2ip is given by 



f{x) = ln„(x) 



(5.23) 



Second functional equation: The following first order 
differential-functional equation emerges in statistical me- 
chanics within the context of the maximum entropy prin- 
ciple 



[x f{x)]=Xf{x/a) 



d_ 

dx 

/(1) = , 



(5.24) 

(5.25) 
(5.26) 
(5.27) 



a and A being two arbitrary constants. The latter 
problem admits two solutions. The first is given by 
f{x) = ln(x) and a = 1/e, A = 1. The second solution 
is obtained after tedious but straightforward calculations 
and is given by 



f{x) = In^ix) , 



and 



l-K 
1 + K 



A = 



. \ l/2ft 



(5.28) 

(5.29) 
(5.30) 



D. The Entropy 

A physically meaningful link between the functions 
\iii^{x) and exp^(x) is given by the following variational 
principle. 

Let be h(q) an arbitrary real function and f{q) a real 
positive function of the variable q E A. The solution of 
the variational equation 



dq f{q) In, f{q) + / dq f{q) h{q) 



is unique and is given by 

f{q) = a cxp^{h{q)/X) , 



= , 
(5.31) 

(5.32) 



a and A being the constants defined by Eqs. (|5.29p and 
(I5.30p . The solution of the variational equation (|5.3ip is 
trivial and employs Eq. (|5.24p . 

This important result permits us to interpret the func- 
tional 



5k = - / dq f{q) In^ f{q) , 
J A 



(5.33) 



which can be written also in the form 

/(g)^' "-/(g)^ +" 
/A 2k 



Sk = dq 
J A 



(5.34) 



as the entropy associated to the function exp^, (x) . It 
is remarkable that in the k ^ limit, as \t\ii{x) and 
exp, (x) approach ln(x) and exp(x) respectively, the new 
entropy reduces to the old Boltzmann-Shannon entropy. 

It is shown that the entropy has the standard prop- 
erties of Boltzmann-Shannon entropy: is thermodynami- 
cally stable, is Lesche stable, obeys the Khinchin axioms 
of continuity, maximality, expandability and generalized 
additivity. 



VI. OTHER REPRESENTATIONS OF THE 
K-EXPONENTIAL 

A. Rapidity variable 

We define the rapidity as r = c arctanh (w/c), so that 
in the classical limit it reduces to the particle velocity. 
The rapidity in dimensionless form p = r/cn^ is given by 



p(u) — — arctanh(Kw) 

K 



(6.1) 



Clearly the above expression of the rapidity is written 
in terms of the velocity, but it can be written also in 
terms of other variables. For instance, after taking into 
account the standard formulas of relativistic dynamics 
e.g. u = qj -^/l -f K^q^ etc, the rapidity can be expressed 
in terms of the momentum, of the total energy, and of 
the kinetic energy, according to 



p{q) = — arcsinh(Kg) 

K 

p{£) — — arccosh(K^£) 

K 

p(W) = - arccosh(l -I- k^W) , 



and obviously it results in 



p = p{q) = p{u) = p[£) - p{W) 



(6.2) 
(6.3) 
(6.4) 

(6.5) 



B. 



The function K-exponential in the velocity 
representation 



We recall that the composition laws of the momenta 
and of the velocities are two generalized sums. For this 
reason we will consider more in detail the rapidity vari- 
able in the momentum representation and in the velocity 
representation. From Eq. (|6.5p we have 

(6.6) 
(6.7) 



exp(p) = exp(p(g)) = exp(p(u)) 
or equivalently 

exp(p) = exp„(g) = exp''(u) . 
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The functions exjp^(q) and exp'^(w) are defined through 



exp^ (q) = exp ( — arcsinh Kg 



(6.8) 



exp''(w) = exp ( — arctanhKii ) , (6.9) 



or equivalently 



expjg) = V1 + «V + 



Kq 



1/k 



exp''(w) 



1 + KU 
I — KU 



1/2k 



(6.10) 
(6.11) 



The explicit relationships linking exp^ (q) , and exp'^ (u) 
are given by 

exp"(u) = exp, (u7(u)) , (6.12) 
expjg) = exp" (g/7(g)) , (6.13) 

7 being the Lorentz factor linking momentum and veloc- 
ity through g = 7 u, which has the following expressions 
in the velocity representation and in the momentum rep- 
resentation 



7(u) 



1 



Vl — K^u'^ 



7(g) = + K^q^ , 



(6.14) 
(6.15) 



respectively. 

We can conclude that the three functions exp(p), 
exp^(g), and exp'^(M) are the same function in three dif- 
ferent representations. 

Hereafter we discuss briefly the function exp''(u) which 
has been introduced firstly in the appendix of ref. Jjl]. 
Starting from the velocity generalized sum we can obtain 
easily the K-differential in the velocity representation as 
follous 



du 



1 — K^U^ 



(6.16) 



The K-derivative and the K-integral of the function f{u) 
with respect the dimensionless velocity, assume respec- 
tively the forms 



d'^u 



= (1 - 



df{u) 



du 



d^u f{u) = 



du 



1 — K^U^ 



fin) 



(6.17) 



(6.18) 



The function f{u) ~ exp'^(it) can be obtained as solu- 
tion of the two following equations 



f{ui ® U2) = f{ui)f{u2) , 



dfju) 
d'^u 



fin) 



(6.19) 



(6.20) 



It is easy to verify that exp''(M) is a monotonic, contin- 
uous function and is defined in the interval — 1/|k| < u < 
1/\k\. The inverse function of exp''(u) namely ln'"(a;) is 
given by 



In'" (a;) = — tanh(Klna;) 

K 

In'" (a;) = — 



(6.21) 
(6.22) 



The relationships linking ln'*'(a::) and hii^{x) are the fol- 
lowing 



ln''(a;) 
ln„(x) 



ln„(x) 



^Jl + ti^\n^{xf 

ln"(a;) 
^l-K^\n^ix)]' 



(6.23) 
(6.24) 



The properties of the functions exp''(x) and In'" (a;) fol- 
lows easily from the ones of the functions exp^(a;) and 
In^lx). For instance, we obtain 



exp''(a;) exp''(— x) = 1 
ln'^(l/x) = -ln'^(a:) , 



(6.25) 
(6.26) 



and so on. 

In ref. [l^ a general procedure to deform the ordinary 
mathematics starting from an arbitrary generalized sum 
has been proposed. This procedure have been adopted in 
ref. [IBl to construct the deformed mathematics related 
to the momentum generalized sum and then, based on 
the function exp^(g). 

Clearly starting from the velocity generalized sum we 
can construct the related deformed mathematics based 
on the function exp''(u). 

For instance we can introduce the deformed hyperbolic 
functions sinh''(u) and cosh'^(w) according to 

3inh"(,) = exp-(.)-^exp-(-u) ^7) 

cosh'^(u) = ^^^P^M±^^^P^ . (6.28) 
2 

These two functions defines a K-deformed hyperbolic 
trigonometry which is isomorphic with respect the or- 
dinary hyperbolic trigonometry. 

The deformed cyclic fimctions sin''(u) and cos''(u) 
given by 

. ^ exp'^(m) -exp'^(-zu) 
sm u = , 6.29) 

^, , exp''(zu) -I- exp''(-zu) ,^ „r.\ 
cos'^(u) = ^ ' ^ ^ , (6.30) 

defines the K-deformed cyclic trigonometry which is iso- 
morphic with respect the ordinary cyclic trigonometry. 

The three mathematical structures based on the ordi- 
nary sum, on the momentum generalized sum, and on 
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the velocity generalized sum, result to be mutually iso- 
morphic. 

Regarding the possibility to consider further represen- 
tations of the exponential function in special relativity 
we recall that the composition laws of the total energy 
£, and of the kinetic energy W, are not generalized sums 
and for these reason we can't use the variables £ and W 
to introduce new representation of the exponential func- 
tion. Let us consider for instance the function 

/(W) = exp ( p{W)) = exp, ( g(W)) . (6.31) 

Clearly this function can't be considered as a represen- 
tation of the exponential function. Indeed it results in 

/(W)/(->V)^l . (6.32) 

In general, a deformation of the ordinary exponential 
i.e. exp {h{x)) is not an its representation. In order to 
have a new representation of the exponential (the mathe- 
matics underlying the new exponential, must be isomor- 
phic with respect the ordinary mathematics, underlying 
the ordinary exponential) the function h{x) must have 
specific properties [18| . 

It is important to note that in the classical limit, the 
composition law for the relativistic kinetic energies, given 
by Eq. (|2.16p . reduces to the following expression of 
Galileian relativity 

Wi © W2 = Wi + W2 + ■ (6.33) 

Then, a non trivial composition law, which is not a gen- 
eralized sum, appears also in classical physics when we 
change the particle observation inertial frame. The latter 
composition law never has been used in classical physics 
to introduce new representations of the exponential func- 
tion. We stress that in classical physics, the ordinary 
exponential can be viewed as a function emerging in the 
momentum or equivalently in the velocity representation. 
Indeed, according to the Galileian relativity the compo- 
sition laws for the classical momenta and velocities are 
the ordinary sum and this ordinary sum generates the 
ordinary exponential and the ordinary mathematics. 

In conclusion in special relativity emerge only two non 
trivial representations of the exponential function namely 
the deformed exponentials exp^(g) and exp''(M) defined 
by Eqs. (|6.10p and (|6.1ip respectivelly. There are sev- 
eral reasons to select the function exp^(g) as the more 
proper representation of the exponential function in or- 
der to construct the relativistic statistical mechanics. 

A first reason is related to the fact that the momentum 
is the more proper variable to formulate the one-particle 
relativistic dynamics. Indeed the dynamical Lorentz 
transformations involve the momentum and no the ve- 
locity. Also the relativistic Newton equation in terms of 
the momentum, assumes a very simple form, similar to 
the classical Newton equation. On the contrary the rela- 
tivistic Newton equation, if expressed in terms of velocity 
and acceleration, results to be a non-linear equation very 
different with respect the classical Newton equation. 



A second reason is related to the fact that also in the 
formulation of the many body relativistic theory the mo- 
mentum results to be the more proper variable. For in- 
stance in relativistic kinetics but also in relativistic field 
theories, the Lorentz invariant integration involves the 
momentum and no the velocity. 

The third and more important reason is related to the 
maximum entropy principle, the cornerstone of statisti- 
cal mechanics. Indeed, the couple of functions exp^((7) 
and \nK,{q) are linked through the variational principle 
described in the subsection D of the section V. On the 
contrary the functions ex.'p'^{q) and ln''(g) can't be con- 
nected by the same variational principle. 

VII. RELATIVISTIC KINETICS 

A. The K-Product and K-Sum of functions 

Let us consider the set of the non negative real func- 
tions V = {/, h, w, ...}. 

Proposition 1: The composition law ® defined through 

K 

ln«(/® /i) = ln,/ + ln,/i , (7.1) 

which reduces to the ordinary product as k ^ 0, namely 
f ® h = f ■ h, is a, generalized product and the algebraic 

structure (T> — {0}, (8)) forms an abelian group. 

Proof: Indeed this K-product has the following proper- 
ties 

1) associative law: (/ ® h) 1^ w — f (E) {h (g) w); 

2) neutral element: f®l — l'E>f = f; 

3) inverse element: / ® (1//) = [1/ f)® f = 1; 

K K 

4) commutative law: f h — h® f . 

Of comse the Ae-division can be defined as follows 
f(Z)h^f®il/h). 

K K, 

The deformed A^-power f^'^ is defined through 

ln„(r'-) =rln„/ , (7.2) 

and generalizes the ordinary power /''. In particular, 
when r is integer one has f^^ = / ^ /••• ^ /, (r times). 

K K 

The K-product allows us to write the following property 
of the K-exponential 

exp^(a:) (g) exp^(?/) = exp^(a; + y) . (7.3) 
Proposition 2: The algebraic structure {V,®) forms 

K 

an abelian monoid. 

Proof: Indeed the element does not admit an inverse 
element. 

Furthermore, just as in the case of the ordinary product, 
it results / (g)0 = 0«)/ = 0. 

Proposition 3: The composition law © defined through 

exp ( ln« (/ e /i) ) = exp ( ln« /) + exp ( ln« h) , (7.4) 
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which reduces to the ordinary sum as the deformation 
parameter approaches to zero, namely / = / + isa 



generahzed sum and the algebraic structure (V, ®) forms 
an abelian monoid. 

Proof: Indeed this K-sum has the following properties 

1) associative law: {f (B h) Q) w = f (B {h® w); 

2) neutral element: / ©0 = 0e/ = /; 

K K 

3) commutative law: / © h — h(B f. 

K K 

Proposition 4: The product (E) and sum © are distribu- 

K K 

tive operations 

w® [f ® h) ^ {w® f)® {w® h) . (7.5) 



B. Evolution Equation 

By using the standard notations of the relativistic the- 
ory we denote with x — x'^ = {ct, x) the four-vector 
position and with p ^ = {p'^,p) the four- vector mo- 
mentum, being — -\- rn?-(? and employ the metric 

Let us consider the following relativistic kinetic equa- 
tion 



we have 



S'p' d^pi d^p[ 



G 



x[C{f'J[)-C{f,h)], (7.6) 

where the distribution / = f{x,p) is a function of the 
four-vectors position and momentum, G is the transition 
rate which depends only on the nature of the two body 
particle interaction and C(/, /i) is the two particle corre- 
lation function with the same four-vector position x, and 
four-momenta p and pi respectively. 

We note that the left hand side of Eq. (|7.6[) is the 
same as in the standard relativistic Boltzmann equation. 
In the particular case where the two particle correlation 
function is assumed to have the same expression like in 
the classical Boltzmann equation i.e. C{f,fi) = f fi 
(Stosszahlansatz), the above equation reduces to the or- 
dinary relativistic Boltzmann equation (2l| . admitting 
as stationary distribution an exponential distribution, 
known as relativistic Maxwell-Boltzmann distribution or 
as Juttner distribution. 

Clearly in the case where C(/,/i) 7^ //i, Eq. 
describes a new relativistic kinetics, radically different 
from the standard one. In the following we pose 



C(/,/i) = (//«) ® {hi a) 



(7.7) 



The origin and value of the constant a will be discussed 
in the following. 



(//«)®(/i/«) = (/'/«) 



(7.1 



and after taking into account the definition (|7.ip of the 
K-product, we obtain 

ln,(//a) + InJ/i/a) = ln,(/7«) + WU'i/a) ■ (7.9) 

This last equation represents a conservation law and then 
we can conclude that ln^{f/a) is a summational invari- 
ant; in the most general case it is a linear combination 
of the microscopic relativistic invariants, namely a con- 
stant and the four-vector momentum. In the literature 
it is shown that in presence of external electromagnetic 
fields the more general microscopic relativistic invari- 
ant has a form proportional to {p^ + qA^/c) U,^+ con- 
stant, being the hydrodynamic four- vector velocity 
with U^Uu = c^. Then we can pose 

In, (//a) ^^^[{p- + qA'^/c) U,~m^]+^. (7.10) 

Consequently we obtain the following stationary distri- 
bution 

/ = aexp, I ~P 1 . (7.11) 

At the moment a. A, /?, and remains arbitrary con- 
stants which will be calculated and/or interpreted by us- 
ing the Maximum Entropy Principle imposing that the 
stationary distribution of the system (|7.1ip must maxi- 
mize the entropy of the system. 



D. The maximum Entropy Principle 

We define the four-vector entropy S" — (5*°, S) as fol- 
lows 



d'^p 



(7.12) 



The identity d^p/p^ = d^p 2 9{p^) 5{p^p^i — m?c^) per- 
mits us to write S'^ also in the form 

S'' = - [ d*p2e{p°)S{p^'Pf,-m^c^) p\f In,/ . (7.13) 



In the latter expression d^p is a scalar because the Jaco- 
bian of the Lorentz transformation is equal to unit. Then 
since transforms as a four- vector, we can conclude that 
S'^ transforms as a four-vector. 

The quantity S is the entropy flow while 5° = 5* is the 
K-entropy given by 



C. Stationary distribution 



S = - d'pflnj 



(7.14) 



We consider now the steady states of Eq. (|7.6p for 
which the collision integral becomes equal to zero. Then 



The maximization of the latter entropy under the con- 
straints imposing the conservation of the norm of the 
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distribution /, and the a priori knowledge of the value of 
the more general microscopic invariant, conducts to the 
following variational equation 



£p f \nj + f3fi / d^p f 



Sf 



-/3 J d'p [ [p" + qA'^/c) - mc'\ / j - , (7.15) 

/3 and ji being the Lagrange multipliers. 

The solution of the latter variational problem conduct 
to the stationary distribution (|7.1ip only thanks to the 
fact that the function IukX/") has the property 



df 



[/ln,(/)]-Aln,(//a) , 



with 



1 - K 

1 + K 



. \ l/2« 



(7.16) 

(7.17) 
(7.18) 



It is important at this point to remark that the func- 
tion In" (a;) does not possess the latter property. For this 
reason, we can't use the generalized logarithm and ex- 
ponential in the velocity representation, in order to con- 
struct a relativistic statistical theory. 

The distribution function (j7.1ip in the global rest 
frame where Ui, = (c, 0, 0, 0) and in absence of external 
forces i.e. — Q simplifies as 



f = a exp^ -13 



W- n 
'J 



(7.19) 



being W the relativistic kinetic energy. 

We observe that the latter distribution in the classical 
limit {k ^ 0, W ^ 0) reduces to the classical Maxwell- 
Boltzmann distribution i.e. / « (1/e) exp (— /3 {W — fJ,)), 
while at relativistic energies {W —>■ +oo) presents power 
law tails / oc W~^^'^, in accordance with the experimen- 
tal evidence in several relativistic systems. 



E. The H-theorem 

In the ordinary relativistic kinetics it is well known 
from the H-theorem that the production of entropy is 
never negative and in equilibrium conditions there is no 
entropy production. In the following we will demonstrate 
the H-theorem for the system governed by the kinetic 
equation (|7.6p when the two-particle correlation function 
is given by (|7.7p . 

By using the property (|7.16p of lnK(/) and the notation 
g = f /a one obtains 



a, (/In,/) = 



■^7 f li^'' / 



XW{f/a)d,f 
Xaliii^{g)d^g . 



(7.20) 



The entropy production di.S'^ can be calculated start- 
ing from the definition of S'^ given by Eq. (|7.12p . By 
using the result (|7.20p and the evolution equation (|7.6p 
it obtains the following expression for the entropy pro- 
duction 



= -A 



Ao^y" ^ In^Cg) p" dug 
d'^p 



p 



I 

^"^nig) p" d^f 



^ d^p' d^pi d^p[ d^p ^ 



p'O pO p'O po 

X (g' (^g'l -g^gi) ^^^ig) 



(7.21) 



Since the Lorentz force F'^ has the properties p^F^ ~ 
and dF'^/dp'^ = the last term in the above equation 
involving F'^ is equal to zero [2l[, therefore we have 



-A 



G 



d^p' d^pi d^p'i d^p 

X (g' ®g'i- g^gij in,.(5)- 



(7.22) 



Given the particular symmetry of the integral in the lat- 
ter equation we can write the entropy production as fol- 
lows 



d^p' d^pi d^p[ d^p 



G 



X (g' (^g[ - g ®5i ) 

X [ln,(.g) + ln«(5i) - ln,(.g') - ln,(.gi)] .(7.23) 
Finally we set the latter equation in the form 



d.S" = - A 



d'^p' d^pi d^p[ d^p 



G 



y Pi 
ig' ®g'i-g® gi 

\ K K 

In, (g' (g)5i) - In, (5(8)51) 



(7.24) 



and after posing h = g ®gi and h' = g' ®g'i and taking 

K K 

into account that In, h is an increasing function, we note 
that {h' - h) (In, h' - In, ft.) > 0, V /i, h' . 
Consequently we can conclude that 



d^S" > 



(7.25) 



This last relation is the local formulation of the relativis- 
tic H-theorem which represents the second law of the 
thermodynamics for the system governed by the evolu- 
tion equation (|7.6p and the two-particle correlation func- 
tion (fTT)) . 
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F. On the definition of tlie entropy 



VIII. CONCLUSIONS 



Let us adopt for the two particle correlation function 
the definition 



(7.26) 



in place of the (|7.7p . It is straightforward to verify that 
one obtains the following expressions for the distribution 
function, the four-vector entropy and the scalar entropy 

f = exp^ - P ^ , (7-27) 



(7.28) 
(7.29) 



S^ - fd^pf ln,(a/) , 



in place of ^LU\ . flAQ^ and (fTM]) respectively 

In particular, in the classical limits we have liniK— ►o — 
1/e, and the entropy l|7.29p reduces to the classical ex- 
pression Scias = — J d^p f [\nf — 1 ] , used some time in 
kinetic theory [2l| . 

In order to explain better the differences between the 
two definitions (|7.14p and (|7.29p of the entropy we con- 
sider a system described by a discrete probability dis- 
tribution / = {/i , 1 < i < N}. We observe that the 
probability distribution / = {Si„i , 1 < i < N} being m 
a fixed integer with 1 < m < iV, describes a state of 
the system for which we have the maximum information. 
It is natural to set for this state S = 0. Clearly it is 
possible only if we adopt for the system entropy the def- 
inition (I7.14p . This interesting property of the entropy 
(|7.14p makes it more appealing with respect the defini- 
tion (|7.29p which predicts a residual entropy 5 = — In^ a 
for the state corresponding to the maximum information. 
Anyway, the two above discussed choices for the entropy 
definition does not influence the physics of the system. 



We have shown that the special relativity laws and the 
maximum entropy principle, suggest a relativistic gener- 
alization for the two-particle correlation function in the 
relativistic Boltzmann equation. This fact imply a rela- 
tivistic generalization of the classical Boltzmann-Gibbs- 
Shannon entropy. 

The so obtained, fully relativistic Boltzmann equation, 
obeys the H-theorem and predicts a stationary stable dis- 
tribution, presenting power-law high-energy tails, accord- 
ing to the experimental evidence. The ensued relativistic 
kinetic theory preserves the main features of the classical 
kinetics which recovers in c ^ oo limit. 

In the last few years the statistical theory based on 
the new entropy [H, [13, d, [Si , has been considered by 
various authors. Investigations related with the founda- 
tions of the theory include e.g. the H-theorem and the 
molecular chaos hypothesis [12, , the thermodynamic 
stability [2l[2i], the Lesche stability [2i,|23,[2l|2i|, the 
Legendre structure of the ensued thermodynamics (soj 
etc. On the other hand, specific applications of the the- 
ory, include e .g. the cosmic rays [ll|, relativistic [3l| 
and classical [32| plasmas in presence of external elec- 
tromagnetic fields, the relaxation in relativistic plasmas 
under wave- par ticle interactions [s^ [s^l, astrophysical 
systems [Sa], the kinetics of interacting atoms and 
photons l37l|, p article systems in external conservative 
force fields '385, the quark-gluon plasma formation [s^ 
etc. Other applications regard dynamical systems at the 
edge of chaos [4l| , fractal systems [i^l , the random 
matrix theory [43] , the error theory [31 , the game theory 
45], the Information theory etc. Also applications 
to economic systems have been considered e.g. to study 
the personal income distribution [i^, [3, to model de- 
terministic heterogeneity in tastes and product differen- 
tiation m, m etc. 
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